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Abstract 



I An impact parameter representation for soft gluon radiation is applied to obtain both the initial 
P-i! 

Q_i' decrease of the total cross-section (atot) for proton-proton collisions as well as the later rise of atot 
(D ■ 

. with energy for both pp and pp. The non-perturbative soft part of the eikonal includes only limited 

- ^ , low energy gluon emission and leads to the initial decrease in the proton-proton cross- section. On 

$H ' the other hand, the rapid rise in the hard, perturbative jet part of the eikonal is tamed into the 

! 

experimentally observed mild increase by soft gluon radiation whose maximum energy rises slowly 
with energy. 
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I. INTRODUCTION 



In this paper, we extend our approach to the role pl^ed by soft gluon radiation in 
determining the energy dependence of total cross-sections 1], by including new low energy 
features. The experimental information on the total cross-section is now sufficient to allow for 
definite progress towards its description through QCD. Thanks to the recent measurements 
at HERA and LEP 0, Q] providing knowledge of total hadronic cross-sections 

involving photons in the energy interval ^/s =1-7- 100 GeV, we now possess a complete 
set of processes to study in detail and in depth, namely pp, pp, ■yp and 77. The purejy 
hadronic processes are well measured over an extended range, up to cosmic ray energies 
leading to quite precise parametrizations [9'], while the other two allow for probes of the 
hadronic content of the photon versus that of the proton. They also allow for checks of the 
Gribov-Froissart factorization hypothesis jiol. 

The three basic features which the data exhibit and which require a theoretical explana- 
tion and understanding are: (i) the normalization of the cross section, (ii) an initial decrease 
and (iii) a subsequent rise with energy. All three are reasonably well described in the Regge 



12 1 as a sum 



trajectory language, which suggests a parametrization of all total cross-sections 
of powers of the square of the cm. energy s. The oldest and simplest of these parametriza- 
tions is 

atot{s)=Xs^ + Ys-^. (1) 

In this model, the initial decrease refiects the disappearance (with increasing energy) of a 
Regge trajectory exchanged in the t-channel, with rj = 1 — ap{0) where the intercept (at 
t = 0) of the leading Regge trajectory is ap{0) ~ 0.5. At the same time, the rise in the 
cross section is attributed to the exchange of a trajectory (the Pomeron) with the quantum 
numbers of the vacuum, such that ap{0) = 1 + e. e is expected to be a small number so as 



not to defy too much the Froissart bound, and, phenomenologically, e = 0.08 -^ 0.12|l3l \14\. 

The two terms in which the cross-sections are split in Eq.(^ describe well the main 
features of the observed total cross-section data, with an apparently constant behaviour 
at ^/s ^ For proton-proton and proton-antiproton scattering, this occurs between 15 
and 25 GeV, where the cross-section is about 40 mb. At large energies, the Regge term 
disappears, while the first term is important everywhere, and dominates asymptotically. 
Thus the decrease is due to the Regge term, the rise described by the Pomeron with some 
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taming due to the Regge term and the normahsation is determined by both the Reg ge a nd 
the Pomeron terms. More recent analyses add further terms such as log'^{s/ sq) terms 15|. 

In the following we shall discuss these features one by one within the context of QCD, 
using the simple but useful parametrization of Eq.(0) as reference. We shall modify it and 
through it develop further the model of ref. p] to increase our understanding of the energy 
dependence within QCD. In section II, we discuss the QCD origin of the two component 
model of Eq. (P). In section III, we summarise how the Born term of the QCD contribution, 
the mini-jets, is embedded in the eikonal formalism which provides a unitarised descrip- 
tion of the total cross-section. In section IV, we discuss the analyticity requirements upon 
the impact parameter amplitudes of the eikonal formalism. In section V, we review the 
soft gluon transverse momentum distribution on which our model for the b-distribution of 
partons is based. In section VI we discuss the main features of the model 1] for the energy 
dependence of the impact parameter distribution, induced by soft gluon emission. In section 
VII, through soft gluon emission, we obtain the transverse parton overlap function for the 
non-perturbative part of the cross-section. In section VIII, we introduce the normalisation 
of the cross-section. We present our model for atot for pp and pp and compare them with 
currently available data. In section IX, we present an interesting feature of our model: in a 
proton or antiproton collision, the (average) distance between the scattering centres (i.e. the 
constituents) in the transverse space decreases as the energy increases. Similar results have 
been previously obtained in an analysis of the hadronic events at the Tevatron, through an 
impact parameter picture. We note again that our model of the energy dependent impact 
parameter distribution offers a reasonable understanding of both the initial fall and later 
rise with energy of the a tot) whereas for the normalization there is still further work to do. 

Before beginning the detailed discussion, we indicate below briefly and qualitatively, why 
we believe that it is both important and necessary to include effects of soft-gluon emission 
in the mini-jet formalism jlfil . IitI lisj ]. We recall here that in the mini-jet picture the rise 
in cross-sections is driven by the increasing number of low-x gluon-gluon collisions and that 
the predicted rate of the rise is generally found to be uncomfortably large. In section VI we 
give the details of the mechanism by which the effect of soft gluon emissions can reduce the 
rate of this rise, at any given cm. energy. Below we summarise the effect qualitatively. 

In a calculation in the QCD improved parton model, the effect of gluon radiation on the 
longitudinal momentum carried by the partons is included, at least in part, in the Dok- 
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shitzer, Gribov, Lipatov, Altarelli and Parisi (DGLAP) (l9| evolved parton distribution. 
Here we crucially look at the effect of initial state gluon radiation on the transverse momen- 
tum distribution of partons. In our model the soft gluon resummed transverse momentum 
distribution of partons in the hadrons and the parton distribution in impact parameter 
space are Fourier Transforms of each other. The larger the transverse momentum the larger 
is the acoUinearity of the two colliding partons, leading to a reduction in parton luminosity 
and hence to a reduction of the cross- sect ion. The higher the cm. energy of the parent 
hadron, the more energetic are the parent valence quarks emitting gluons and the more is 
the acoUinearity of the two partons involved in the parton subprocess. Hence, we expect the 
effect to be dependent on the cm. energy. 



II. THE TWO COMPONENT MODEL 



Before entering into the technical details about unitarization, it is well to ask (i) where 
the "two component" structure of Eq.(^ comes from and (ii) why the difference in the two 
powers (in s) is approximately a half. 

Let us first remember that the two terms of Eq. (Q) reflect the well known duality between 
resonance and Regge pole exchange on the one hand and background and Pomeron exchange 
on the other, established in the late 60's through FESR [20]. This correspondance meant 
that, while at low energy the cross-section could be written as due to a background term 
and a sum of resonances, at higher energy it could be written as a sum of Regge trajectory 
exchanges and a Pomeron exchange. 

Our present knowledge of QCD description of hadronic phenomena gives further insight 
into the nature of these two terms. We shall start answering the above two questions 
through considerations about the bound state nature of hadrons which necessarily transcends 
perturbative QCD. For hadrons made of light quarks(g) and glue(5'), the two terms arise 
from qq and gg excitations. For these, the energy is given by a sum of three terms: (i) the 
rotational energy, (ii) the Coulomb energy and (iii) the "confining" e nerg y. If we accept 
the Wilson area conjecture in QCD, (iii) reduces to the linear potential 2l|, |22|. Explicitly, 
in the CM frame of two massless particles, either a gg or a gg pair separated by a relative 
distance r with relative angular momentum J, the energy is given by 

HJ,r) = ^-^ + arr, (2) 



where i = 1 refers to qq, i = 2 refers to gg, t is the "string tension" and the Casimir's 
are Ci = Cp = 4/3, C2 = Cq = 3. a is the QCD couphng constant evaluated at a 
some average value of r and whose value will disappear in the ratio to be considered. The 
hadronic rest mass for a state of angular momentum J is then computed through minimising 
the above energy 



M,(J) = Minr[^ - ^ + C^rrl (3) 



which gives 



Mi{J) = 2^ {Qt)[2J - Qa]. (4) 

The result may then be inverted to obtain the two sets of linear Regge trajectories ai{s) 

C ^ X 

ai{s) = ^ + {-^^)s = ai{0) + a'iS. (5) 
Thus, the ratio of the intercepts is given by 

Using our present understanding that resonances are qq bound states while the background, 
dual to the Pomeron, is provided by gluon-gluon exchangesj^], the above equation can be 
rewritten as 

If we restrict our attention to the leading Regge trajectory, namely the degenerate p — uj — cf) 
trajectory, then q;_r(0) = rj ^ 0.48 — 0.5, and we obtain for e ~ 0.08 — 0.12, a rather 
satisfactory value. The same argument for the slopes gives 

^ = C,/Ca = i (8) 

SO that if we take for the Regge slope a'p f» 0.88 — 0.90, we get for a'p 0.39 — 0.40, in 
fair agreement with lattice estimates 24| . 

We now have good reasons for a break up of the amplitude into two components. To 
proceed further, it is necessary to realize that precisely because massless hadrons do not 
exist, Eq.(PQ) violates the Froissart bound and thus must be unitarized. To begin this task, 
let us first rewrite Eq.((T)) by putting in the "correct" dimensions 

atot{s) = a,{s/sy + a2is/sy/^, (9) 
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where we have imposed the nominal value t] = 1/2. In the following, we shall obtain rough 
estimates for the size of the parameters in Eq.Q. 

A minimum occurs in atotis) at s = s, for (T2 = 2eai. If we make this choice, then 
Eq.Q has one less parameter and it reduces to 

atotis) =a,[{s/sY + 2e{s/sY/']. (10) 

We can isolate the rising part of the cross-section by rewriting the above as 

atotis) =ai[l + 2eis/s)'/^] + a,[is/s)' - 1]. (11) 

Eq.llllj) separates cleanly the cross-section into two parts: the first part is a "soft" piece 
which shows a saturation to a constant value (but which contains no rise) and the second a 
"hard" piece which has all the rise. Morover, s naturally provides the scale beyond which 
the cross-sections would begin to rise. Thus, our "Born" term assumes the generic form 

^tot("S) = (^softis) + 'dis - s)ahardis). (12) 

with asoft containing a constant ( the "old" Pomeron with ap(0) = 1) plus a (Regge) 
term decreasing as l/^/s and with an estimate for their relative magnitudes (o"2/o"i ~ 2e). 
We shall assume that the rising part of the cross-section atard is provided by jets which are 
calculable by perturbative QCD, obviating (atleast in principle) the need of an arbitrary 
parameter e. 

An estimate of ai may also be obtained through the hadronic string picture. Eq.(j21) 
gives us the mean distance between quarks or the "size" of a hadronic excitation of angular 
momentum J in terms of the string tension 

fiJ)' = ^1^1^, (13) 

r 

So the size Ri of the lowest hadron (which in this Regge string picture has J = 1, since 
aniO) = 1/2) is given by 

Rl = - = 8a' (14) 

T 

If two hadrons each of size Ri collide, their effective radius for scattering would be given by 



Reff = jRl + Rl = V2R,, (15) 
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and the constant cross-section may be estimated (semi-classically) to be roughly 

ai = 2nRlff = AttRI ~ — = 327ra', (16) 

which is about 40 mb, a reasonable value. In the later sections, for the "soft" cross-section 
we shall take a value of this order of magnitude as the nominal value. 

The unitarization now proceeds very simply by eikonal exponentiation in impact param- 
eter space, as described in detail in the sections to follow. 



III. THE RISE WITH ENERGY 



QCD offers an elegant explanation o 



out by several authors in the past 
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the rise in the minijet formalism as has been pointed 
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18|. The original suggestion was that the rise of 



atot with energy is driven by the rapid rise with energy of the inclusive jet cross-section 

dpt I dXi [ dX2 fi\a{Xl)fj\b{x2) ^^'^~^^^^^\ (17) 



where subscripts a and b denote particles (7, p, . . .), i, j, k, I are partons and Xi,X2 the 
fractions of the parent particle momentum carried by the parton. s = X1X2S and a are hard 
partonic scattering cross-sections. Note that da/dpt oc p^^; the cross-section defined in 
Eq. lfTTj) therefore depends very sensitively on ptmin, which is supposed to parametrize the 
transition from perturbative to nonperturbative QCD. The rise of the inclusive jet cross- 
section with energy is understood in terms of the increasing number of hard partons which 
gives rise to an increasing probability for the occurrence of hard scattering processes. Quan- 
titatively, factorisation of QCD allows us to use the currently available parametrisations 
of the scale dependent parton densities and calculate the energy dependence of the result- 
ing jet cross-section, by convoluting the parton densities with the subprocess cross-section 
determined by perturbative QCD. 

For any fixed ptmin, typically 1^2 GeV, one finds that this cross-section is a steeply rising 
function of energy. If y/s ^ Ptmin, the integral in Eq. lfTTjl receives its dominant contribution 
from xi^2 <^ 1- The relevant parton densities can then be approximated by a simple power 
law, / oc In case of pp or pp scattering, a = b and the cross-section asymptotically 

scales like j^] 

1 I s \ s . . 



if J > 1. For J ~ 1.3, as measured by HERA, the jet cross-section will therefore grow 
much faster than the total pp cross-section, which only grows oc log^ s (Froissart bound 



26[), or, phenomenologically II4] for < 2 TeV, oc Eventually the jet cross- 



section ()17|1 will therefore exceed the total pp cross-section. In fact, this rise is far more 



22|. This 



violent than the experimentally observed gentle rise of the total cross-section 
has led to various phenomenological strategies (among them the eikonal formalism |28j |) 
directed towards softening this rise. The apparent paradox is solved by the observation that, 
by definition, inclusive cross-sections include a multiplicity factor. Since a hard partonic 
scattering always produces a pair of (mini-)jets, we can write 

< = (^jet pairXnel, (19) 

where (njet pair) is the average number of (mini-)jet pairs per inelastic collision, a^^* > a^^^i 
then implies (njet pair) > 1? which means that, on average, each inelastic event contains more 
than one hard partonic scatter. The simplest possible assumption about these multiple 
partonic interactions is that they occur completely independently of each other, in which 
case rijet pair obeys a Poisson distribution. At a slightly higher level of sophistication, one 
assumes these interactions to be independent only at fixed impact parameter 6; indeed, it 
seems natural to assume that events with small h usually have larger njet pair- This leads to 
the eikonal formalism mentioned above. 

Convenient and elegant, the eikonal method reduces the rise of this cross-section, and 
allows one to enforce the requirement of s-channel unitarity. Here one obtains the total 
cross-section through the eikonal formula 

atot = 2 y d%l - e-^^'^^^^'^coslkexih, s)] (20) 

and one introduces the jet cross-section as the term which drives the rise in the eikonal 
function. This can be done unambiguously by defining the inelastic cross-section j^l given 
in the eikonal formulation by 

ainel= j d%l - (21) 

This expression can also be obtained upon summing multiple collisions which are Poisson 
distributed with an average number n{b,s) = 2 53m x(6, s). Making the approximation 
3?ex = 0, one obtains a very simple expression 

atot = 2 / d%l - e-"('''^)/2] (22) 



To proceed further, one needs to introduce the soft processes, which cannot be described by 
perturbative QCD. Following the separation shown in Sect. II, one can approximate n[b, s) 
by introducing a separation between soft and hard processes as: 

n{b, s) = Usoft + rihard = Asoft{b)asoft{s) + Ajet{b)ajet{s) (23) 

The separation between hard and soft processes is of course approximate and so is the 
factorization into energy and transverse dimension dependence. However, it is useful as 
it allows one to break down the calculation into building blocks, which can be separately 
understood and put together again later as part of the whole structure. 

The procedure of Eq.(|^ reduces the fast rise due to the mini-jet cross-section, but the 
extent to which this softens the riseMs highly dependent on the impact parameter (6)- 
dependence of n{b,s). The simplest |23] ansatz, which introduces a minimum number of 
parameters is to assume that the 6-dependence is the same for both the soft and the jet 
component, and further that it is given by the Fourier transform of the electromagnetic 
form factor of the colliding hadron, J^(g')|2^]. Thus for protons and antiprotons one will 
have 

As^b) = AUh) = ^ / <fy-^[:FM? = / (24) 

It is well known that these Eikonal Minijet Models (EMM) are unable to properly 
reproduce the experimentally observed, complete rise of the cross-section from the beginning 
up to the asymptotia, without introducing further parameters. As an example, results 
obtained using Eq.(j2H) with u = 0.71 GeV^, and current Gliick, Reya and Vogt (GRV) 
parton densities for the proton 31] to calculate the jet cross-sections in Eq. ()17|l . are shown 
in FIG. dH). One can see from the figure that a ptmin ~ 2 GeV is needed, in order to obtain 
a numerical value of the total proton cross-sections in the 80 mb range, at Tevatron energies. 
However, for such ptmin = 2 GeV, the rise does not begin until -y/i is already in the 100 GeV 
region. On the other hand, a smaller value for the regulator ptmin, typically just above 1 
GeV, would allow for the beginning of the rise around 20^30 GeV, as the data indicate, but 
then the cross-section rises too rapidly in comparison with the Tevatron data. In FIG. (PJ, 
in all jet cross-sections computed using Eq. (|T7jl . we use the strong coupling constant as at 
scale pt- Even at pt = ptmin, this value for the scale, albeit low, is still in a range where the 
asymptotic freedom expression is approximately valid. For the low energy behaviour, all the 
curves shown in FIG. ^are obtained with the same phenomenological fit as in ref.U^. 
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FIG. 1: Comparison between data |a, |3j, |33, |3J, |35|, |3g, |37| and the EMM (see text), for different 
minimum jet transverse momentum. 



Different solutions to the above problems have been adopted, including energy dependence 

Q, Q with 



or using Eq. 



of the cut-off parameter ptmin, adding more terms in Eq. 
different constants for the low and the high energy part. 

We have taken the approach to keep fixed ptmin (for a given beam and target combination) 
and have an energy dependent overlap function, the energy dependence being modelled by 
a QCD motivated calculation. The QCD motivated model gives a transverse momentum 
distribution of the partons which is energy dependent and this, in turn, makes our overlap 
functions energy dependent. We shall see in the next sections that in this model soft gluon 
emission has an effect very similar to that of an energy dependent ptmin- 



IV. ANALYTICITY REQUIREMENTS ON THE IMPACT PARAMETER DIS- 
TRIBUTION 

Within the context of QCD, the problem described in the previous section, sometimes 
referred to as the soft and hard Pomeron problem, can have (at least) two different origins. If 
complete factorisation between the dependence on the energy ^/s and the impact parameter 



10 



b holds, then the abovementioned soft and hard pomeron prolem may simply be taken to 
be indicative that the s-dependence generated through gluon densities, is wrong. On the 
other hand, given the fact that gluon densities are measured in Deep Inelastic Scattering 
experiments, one may think of an alternative explanation of the probelm. This second, com- 
plementary explanation is that not all the s-dependence is due to the jet-cross-sections and 
there is further energy dependence in the impact parameter distribution. Some general ana- 
lyticity arguments can be invoked to shed light on the above and thus limit the arbitrariness 
in the choice of the function s). 

Consider the elastic scattering amplitude Tei{s,t) (with s and t the usual Mandelstam 
variables) normalized so that the total cross section is given by 

atot = 2 J d%l-e-''^^^''''^cos^ex{b,s)] = (^) '^m T,i{s,t = 0), (25) 



and the elastic differential cross section by 



Tell' (26) 



Consistently with Eq. ()25j) . one then has 

Tei{s, t)=isj d^be''^-^[l - e'^(^'^)]. (27) 
For complete absorption, i.e. 3?e x = 0, we get 

Tei{s, t) = is j d'he"^-\l - e-"('''^)/2] (28) 

Let us briefly examine restrictions imposed on the large b-behaviour of the function n(6, s) 
by the requirements of analyticity. Consider the Fourier transform of the elastic scattering 
amplitude, 

A{s,h) = ^ C dt Tei{s,t)Mb^/^) (29) 
where t = —(f, with q the transverse momentum variable. Equivalently, we have 

Tei{s,t) = ms db''Jo{bV^)A{s,b) (30) 



The finite range of hadronic interactions implies that the partial wave expansion converges 
beyond the physical region, i.e., throughout the Lehmann ellipse. This requires that Tei{s, t) 
be analytic in t up to t = where /i is the pion mass. For positive t, we continue the 
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above expression for = ^W^? with W real and positive. In this "unphysical" region, we 
have 

T^i{s, W'^) = iris J db'^Io{bW)A{s, 6^) (31) 
For large b, Io{bW) ~ , so that for the integral to converge, one needs 



2nbW' 

\A{s,b'^)\ < e-^^" with Wo^2^i (32) 

In the purely absorptive model (3fJe x = 0) thus, 

1 - e""(^'^)/2 < ^-bw./2^ (33) 

and we see that n(6, s) must be bounded at least by an exponential. Stronger (but model 
dependent) constraints arise provided one imposes that the elastic differential cross-section 
exhibit a "diffraction peak" . That is 

T,Ks,t)~/(s)e^W* (34) 

where b{s) is the so-called width of the diffraction peak, which has an observed ( approxi- 
mately) logarithmic s dependence. Then, Eq.()29|) gives 



A{s,b)c^'-^ j dq'Ubq)e- 



2718 



2b{s) 



h2 



e (35) 



which requires a Gaussian fall-off of the amplitude in the impact parameter b, with its scale 
determined by the width of the diffraction peak. In the Regge pole description, 

b{s) ~ a'ln{s/so), f{s) ~ -t(3{s/soy+' (36) 

and 

In the EMM, where the impact parameter distribution is given by the Fourier transform 
of the proton form factor, a model we refer to as the Form Factor (FF) model, one has 

n{b, s) = -— {vb) K^{vb)[asoft + cTjet] (38) 
yoTT 

The modified Bessel functions of the third kind K^{z) are bounded by an exponential at 
large values of the argument, i.e. Kf^{z) ~ \[^^^^{'^ + ^{\)}- We see then that the EMM 
in the FF formulation does satisfy the requirements of analyticity in the Lehmann ellipse. 
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But in the EMM, the observed shrinking of the diffraction peak, corresponding to an energy 
dependent near gaussian fall off, is not present. Even if one were to introduce an ad hoc 
energy dependence (as is often the practice) instead of the constant scale parameter z/, as 
in the FF model, still one would be nowhere near the stronger Gaussian decrease at large 
impact parameter values. This is one reason why the FF model in the eikonal formulation, 
where jet cross-sections drive the rise, fails to provide an adequate description of the overall 
energy dependence of the total cross section, without introducing an ad hoc modification of 
the scale parameters. 



V. THE SOFT GLUON TRANSVERSE MOMENTUM DISTRIBUTION 

The model for the total cross-section presented in this paper is based on the ansatz 
that QCD provides the main processes at work leading to the observed energy rise of all 
measured total cross-sections. An important motivation of this model is to make quantitative 
calculations based on current QCD phenomenology, namely current parton densities with 
their energy momentum dependence, running behaviour of the coupling constant and soft 
gluon resummation techniques. Of all these at present soft gluon resummation is the one 
which presents the toughest technical challenge. 

Let us begin by considering the well known function which describes soft gluon emission 
from a parton-parton pair, namely the soft gluon transverse momentum distribution [ssl. Isojl 

''^^^^^ n(ir,)=/-^e---'^(^) (39) 



(i2Kx ^ J (27r) 

with 

h{h) = J d\{k)[l - e-''-^] =J^ Y: \f'\k)j,,{k)\[l - e-^^--^] (40) 



''O i,j=colors 

73. 



where d ng{k) is the distribution for single gluon emission in a scattering process, and 
the QCD current responsible for emission. The above expression has been wide^ 
study the initial state transverse momentum distribution in Drell-Yan processes 



well as W-production 



42 



43|. 



V used 

mm 



used to 
as 



In the limit of large k±h, one can neglect the exponential term, and the above expression 



reduces to the well known Sudakov form factor 



44l |. namely 

d'^k 

" colors 



h{h) ^ Sib) = I d\ik) = E \fmj,Ak)\ (41) 
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Introducing the running coupling constant and its asymptotic freedom expression, the inte- 
gration from 1/6 up to an upper hmit Q gives [45! 

^^^^ (11 -2n;/3) A2'^Zn(l/62A2) ^^^^ 
For small momenta, however, the above expression is not sufficient to reproduce the 
observed tranverse momentum distribution in various hadronic processes. Thus, an intrinsic 
transverse momentum (of Drell-Yan pairs or W-boson or partons, depending on the physical 
p..oces. „„de. co„.de.at,o.) has .o be ...reduced. TKe function „ow .ead. U 

m = fe'pLnt + S{b) (43) 

where the intrinsic transverse momentum p±int is a constant, of the order of a few 100 MeV, 
parametrized according to the process under consideration. 

In our model, the function in Eq. ()40|) . in addition to describing various hadronic trans- 
verse momentum effects, also plays a major role in total cross-section calculations. Our 
model for the hadronic transverse momentum distributions due to soft gluon resummation 
differs in two major points from what we have just described. Basically, we focus on the 
lower and higher limits of integration in Eq. (j40j) . At the lower limit, since this expression 
refers to soft gluons, we suggest that the correct use of this equation requires to integrate 
the gluon momentum down to /c^ = and avoid the introduction of ad hoc quantities such 
as the intrinsic transverse momentum. At the upper limit, one needs to specify, for each 
given process, how the maximum soft gluon energy is defined. However, in order to use the 
above expression for a believable calculation, a number of points need to be clarified, namely 

(i) whether it makes sense to use a parton picture when the gluon momenta become close 

to zero, with the related question of what is the behaviour of the strong coupling 
constant when one integrates the gluon momentum down to zero 

(ii) whether the emitting particles are quarks or gluons 

(iii) what are the constraints from kinematics upon the maximum gluon momentum 

We shall discuss some of these issues in detail in the sections to come, here we comment 
briefly on these three points. 

Concerning the parton picture, while we use it for the mini-jet contribution to the cross- 
section, soft gluon Bloch-Nordsieck resummation factors out of the LO basic scattering 
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process, and is thus independent of parton densities, involving only the momenta and the 
QCD coupling between soft gluons and the emitting partons. On the other hand when 
k± ^ 0, this coupling is not an observable quantity, since it refers to a single soft gluon 
emission, and one soft gluon is not an observable quantity (only its integrated spectrum is). 
This has two effects: firstly, one needs to use a non-perturbative expression for the QCD 
coupling constant, since the momenta are so small, and secondly only the integral of moments 
of as will matter. The infrared behaviour of as is a matter of speculation. We propose our 
own model, whose justification rests upon a Regge description and on the Richardson type 
potential for quarkonium. As will be described in the next sections, a specific form for as is 
chosen, singular in the infrared limit, but integrable [31 . 

Another issue to address relates to the effect of emission from quarks, valence and sea, and 
from gluons in the parton processes. In this paper, we only deal with emission from the initial 
valence quarks, and use the relevant kinematics with their averages. This approximation is 
justified by the fact that, relative to the emission from the initial valence quarks, emission 
from the gluon legs is to be considered as emission from internal legs, thus subleading in 
infrared terms. A complete calculation should of course include also emission from partons 
other than the valence quarks and hence mostly the low-x gluons. We expect this inclusion 
may eventually increase the softening effect. 

To describe the softening effect quantitatively, one needs to focus on the maximum trans- 
verse momentum allowed to single gluon emission. This quantity is energy dependent, as 
one can easily see using the kinematics of single gluon emission in parton-parton scattering 
of initial cm. energy \/§. For the process 

Pi + Pj ~^ gluon + X (44) 

where s = {pi + and X a final state of given momentum Q, the maximum transverse 
momentum of the gluon is given by 

qmaM = ^{^-^) (45) 

If we consider the state X to be the final (mini) jet-jet system in the inelastic collision 
contributing to the cross-section, then the above expression depends on the parton sub- 
energies and on the final state momentum of the jet-jet system, characterized by a transverse 
momentum pt > Ptmin, where Ptmin is a scale chosen to separate hard and soft processes. In 
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principle, for each subprocess of given s and one should evaluate the function h{b) with 
the above Qmax- In practice, we use a value for the maximum transverse momentum allowed 
to single gluons, which is averaged over the initial and final parton momenta. This is shown 
more explicitely in and in section VI, but the result is that for a given ptmin cut-off in the 
minijet distribution, for the valence quarks, the scale Qmax increases with the cm. energy 
of the hadron-hadron system. This can be qualitatively understood by considering that the 
valence quarks will on the average carry a larger energy and can then shed more soft gluons. 
Thus, in the picture we present for hadron-hadron collisions, as the overal energy increases, 
we have more parton-parton collisions for the same ptmin (since the number of low x-gluons 
increases) but also more energy available to soft gluons both from initial valence quarks 
and from all the hard partons in general, and thus more of a reduction. It is the balancing 
of these two effects which we believe to be responsible for the observed softer rise of total 
cross-sections. 



VI. SOFT GLUON EMISSION AND ENERGY DEPENDENCE IN THE IMPACT 
PARAMETER DISTRIBUTION 

In our previous work, we have advocated that a cure for the difficulty in obtaining the 
early dramatic rise and the softer asymptotic behaviour simultaneously, lies within QCD 
itself; viz., the ubiquitous soft gluon emission accompanying all QCD scatterings which can 
slow down any abrupt rise in the cross-section. To make this quantitative, we put forward 
a model for the impact parameter distribution of partons in the hadrons, based on the 
Fourier transform of the transverse momentum distribution of the soft gluons emitted in the 
collisions, as described in the Block-Nordsieck (BN) summation procedure. This distribution 
is energy dependent simply because the maximum energy allowed to each single emitted soft 
gluon, in turn depends on the energy of the colliding partons. In detail, we have a picture 
of parton-parton collisions at all admissible subenergy values and with a given transverse 
momentum due to initial state radiation. In our model the soft gluon resummed transverse 
momentum distribution of partons in the hadrons and the parton distribution in impact 
parameter space are Fourier Transforms of each other. In principle, this formalism could 
be used to obtain impact parameter dependent parton densities, but our aim in the present 
paper is to obtain a prediction for the total cross-section based on currently used QCD 
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functions and parameters such as parton densities and Aqcd- We thus follow our previous 
proposal Q, 0] to average out the behaviour of partons in their transverse momentum 
variable and arrive at the following expression 



2Qm X{b, S) = n{b, S; qmax,Ptmin) = n^oft + nj^t = Asoft{b)asoft + ^B7v(&, qmax)(^jet (46) 

As mentioned earlier, the eikonal formulation provides a natural framework, in which dif- 
ferent contributions to the total cross-section can be resolved into their various structural 
elements: the rise is incorporated in rijet, and the decrease and normalization in risoft- In 
our previous work, we had parametrized phenomeno logically the soft part, and used per- 
turbative QCD for the jet part. In this paper, we study whether soft gluon summation can 
describe the (experimentally observed) initial decrease in proton-proton scattering. 
For this purpose, we write 

n{b, s) = A'^ljasoft + ^BAfCTiet (47) 

with 

g-fc(fe,s) 



Abn = ^ (48) 



where from Eq. ()4Up we have 



h{b, 5) = TT / ir"-(^ ) V ,J [1 - Mkb)] (49) 

Hmax y Hmax 



and qmax depends on energy and the kinematics of the process j4l|. From Eq. (P3|) . the 
following average expression for qmax was proposed in our previous paper 

_< >- ^ ^ ^ ^^^^^^^^^ (60) 

with Zmin = ^Ptminl {^^1^2) and fi/a the valence quark densities used in the jet cross-section 
calculation. 

M establishes the scale which, on the average, regulates soft gluon emission in the colli- 
sions, whereas Ptmin provides the scale which characterizes the onset of hard parton-parton 
scattering. For any parton parton subprocess characterized by a ptmin ~ 1 ^ 2 GeV , M has 
a loga.ithn,ic _ at .ea=o„ably ,ow ene.gy and an al,no=t constant behavtont at Mgh 
energy The eikonal formalism which we use to describe the total cross-section, incorpo- 
rates multiple parton parton collisions, accompanied by soft gluon emission from the initial 
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valence quarks, to leading order. Notice that in this model, we consider emissions only from 
the external quark legs. In the impulse approximation on which the parton model itself is 
based, the valence quarks are free, external particles. In this picture, emission of soft gluons 
from the gluons involved in the hard scattering, is non leading. As the energy increases, more 
and more hard gluons are emitted but there is also a larger and larger probability of soft 
gluon emission : the overall effect is a rise of the cross-section, tempered by the soft emission, 
i.e. the violent mini-jet rise due to semi-hard gluon gluon collisions is tamed by soft gluons. 
Crucial in this model, are the scale and the behaviour of the strong coupling constant which 
is present in the integral over the soft gluon spectrum. While in the jet cross-section never 
plunges into the infrared region, as the scattering partons are by construction semi-hard, in 
the soft gluon spectrum the opposite is true and a regularization is mandatory. We notice 
however that here, as in other problems of soft hadron physics[48], what matters most is not 
the value of as{0), but rather its integral. Thus, all that we need to demand, is that be 
integrable, even if singular Q] . We employ the same phenomenological expression for as as 
used in our previous works, namely 

"'^ ^^ ~ (33-2iV,)ln[l+p(^)^.] ^''^ 



Through the above, we were able to reproduce the effect of the phenomenologically intro- 
duced intrinsic transverse momentum of hadrons 0|, and more recently obtained a very 
good description of the entire region where the total cross-section rises This expression 
for as coincides with the usual one- loop expression for large values of kj_, while going to 
a singular limit for small k_i_. For p = 1 this expression corresponds to the Richardson 
potential l50| used in bound state problems. We see from Eq.f^^ that p = 1, leads to a 
divergent integral, and thus cannot be used. Notice that, presently, in the expression for 
h{b, s), the masses of the emitting particles are put to zero as is usual in perturbative QCD. 
Thus, for _a convergent integral, one requires p < 1 and the successful phenomenology indi- 
cated in |l| gave p = 3/4. However, more study is needed, especially in the full utilization 
of the Bloch-Nordsieck description, before one can completely define an expression for 
in the infrared limit. A different possibility is to use a so-called frozen as model, for which 
cts{kY) = 27in[(fc^ +a2A2)/A2)] ' Thcsc two cxprcssious lead to very different large-b behaviour 
of the function n{b, s) and, in light of the above discussion concerning the shrinking of the 
diffraction peak, give quite a different s-dependence in the rising region of the total proton 
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cross-section. 

In 0, we presented analytic approximations for the function h[b,s), obtaining, in the 
frozen a. case 



hm h{b, M, A) 

6— >oo 



while for the singular case 



lim h{b, M, A) 

fe— >oo 



2cF<ys 



TT 



1 



ln(2M6) + 2 ln(M6) ln(aA6) - In^ (aAb) 



(52) 



2cFb 



TT 



1(1 -p) 



2 ln(2M6) + 



1 — p 



+ 



1 21„(A«)-1 



where M = q„ 



dmax, «s 271na2' " 33~2Nf- 

imately Gaussian limit results in the singular ag case, but not in the frozen case. Indeed, 
for the singular case, one has, to the lowest order, limfe_»oo n{b, s) ~ e~^^'' i.e. an exact 
Gaussian limit for the Richardson potential, which corresponds to p = 1. This provides the 
theoretical reason why the entire region where the total cross-section rises is well described 
by perturbative QCD (jet cross-section) combined with Bloch-Nordsieck summation with a 
singular as in the infrared region. In contrast, neither the frozen ag model nor, the Form 
Factor (FF) model are successful there. 

Other models for the behaviour of in the infrared region, and studies of the range of 
variability of the parameters used in Eq (f5T| will be presented in a forthcoming publication. 



127r 



127r 



From the above, one can see that an approx- 



VII. THE DECREASE PRIOR TO THE ONSET OF MINI-JETS 

Let us now address the question of the lower energy range, prior to the rise, using the 
same phenomenological and theoretical tools of but abandoning the FF and the frozen 
as models which appear to be inadequate. To study the low energy region, we apply our 
procedure to proton-proton scattering, where the absence of resonances in the s-channel and 
leading Regge poles in the t-channel make the picture remarkably simple. At low energies, 
say before 10 GeV in the proton-proton cm. system, one observes a very soft decrease, which 
converts in a rise at an energy of ~ 15 GeV in the cm. In this low energy region, we know 
that gluon densities are still very small, and that (almost) all hard parton-parton scattering 
takes place among valence quarks : a jet, as defined through Eq.lfTTj). is a few thousands 
of the observed atot- In the break up of x into a soft and a hard part, the parameter 
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Ptmin separates hard and soft processes, namely for pf"'^ °" > ptmin one counts parton-parton 
processes as part of the jet cross- sect ion, whereas for < p^.^^^ the process can be 

counted as part of a soft- Thus, in this region, we can study the contribution of valence 
quark scattering without complications from inelastic gluon-gluon collisions. 

This region then exhibits the effect of soft gluon emission accompanying gluon exchanges 
among the valence quarks. At higher energies, these soft interactions still take place and 
be a substantial part of the cross-section, but they will be shielded by the more dramatic 
behaviour of the perturbative QCD processes, since as the energy increases, smaller and 
smaller x-values of the gluon densities are probed and gluon exchanges among gluons start 
becoming important. Thus we must build a piece of the total cross-section which will survive 
at high energies, but which does not contribute to the rise. To begin with, we start with 
a very simple ansatz : that for proton-proton the cross-section agoft is a constant and the 
slight decrease comes from the straggling, acollinearity effect of soft gluon emission. We 
therefore propose, in first instance, the following expression for the average number of soft 
collisions 

nsoftib, s) = A^°^Vo (53) 

with A^jY calculated through Eqs. ()48p49j) and investigate whether it is possible to find a 
constant o"o and a set of parameters {qmax) which can describe pp scattering at low energy. 
For the soft part, the scale qmax corresponds to the maximum energy allowed to soft gluons 
accompanying scattering with a final parton transverse momentum smaller than ptmin- We 
are dealing with soft emission (for hard gluons pt > Ptmin) and thus we expect qmax not to 
be larger than 10 — 20 % of ptmin- This provides an upper bound for qmax for soft processes. 

The observation is then that for processes contributing to n^o/t, a soft gluon will always 
carry away less energy than for those contributing to Uhard- The question is how much lower 
is the allowed energy. We have proceeded phenomenologically, and found a set of values 
which, as will be shown in the last section, can give an acceptable description for cjpp before 
the rise. These values are shown in Table HI Notice that, in order to reduce the number 
of free parameters, we assume that, at low energies, there is only one value of qmax, for 
both hard and soft processes. However, as the energy increases, the scale characterizing 
the soft processes does not grow as the one for the hard case, the latter obtained through 
the kinematics of jet production in a hard parton parton scattering. Here we should expect 
the scale to start as a slowly increasing function of s, but to become a constant as soon as 
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TABLE I: Average Qmax values used for the impact parameter distribution of the soft part of the 
eikonal 







5. 


0.19 


6. 


0.21 


7. 


0.22 


8. 


0.23 


9. 


0.235 


10. 


0.24 


50. 


0.24 


100 


0.24 



hard processes become substantial for y/s > 10 GeV. This is necessary, i.e. qmax does not 
increase indefinitely, because as the energy available to soft gluons increases, at a certain 
point the soft gluons will become hard and then undergo scattering among themselves. 

Clearly, a soft scale not larger than 240 MeV is consistent with our understanding of 
how a proton is structured, if we attribute asoft, the soft component, as the cross section 
when the scattering protons (and antiprotons) manifest themselves as a quark and a spin 
zero diquark. The point is that the slopes for meson and baryon Regge trajectories are 
justifiably equal only if a baryon is pictured as a quark/diquark system similar to a meson 
as a quark/antiquark system (thus making the string tensions equal). This mode for the 
nucleon is soft and diffusely spread over about a fermij^. Consequently, the soft- gluon 
radiation distribution must be limited (lest it break the system). Thus, we estimate for the 
soft process a qmax ~ 1/(1 fermi) = 0.2 GeV. 

With the value of ptmin which gave a smooth description of the total cross-section in Q|, 
we plot in FIG. (0) the behaviour of Qmax as a function of energy, where the upper curve 
is the one obtained using Eq. (j3n|) . for ptmin = 1-15 GeV, whereas the soft qmax starts with 
the same value obtained for the jet term, and then is made to become a constant when it 
reaches 240 MeV. 

With these values for qm.ax we can now calculate Abn for both the hard and soft terms 
in the eikonal as a function of the impact parameter b. Both soft and hard Abn are shown 
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FIG. 2: The maximum transverse momentum allowed (on the average) for single soft gluon emission 
as a function of the cm. energy of scattering hadrons. 

in Fig. 01 for a set of representative cm. energies, y/s = 5, 10, 50, 100, 500 GeV. Notice that 
y4gj^* does not change for y/s > 10 GeV, since Qmax remains constant. 

VIII. NORMALIZATION AND TOTAL CROSS-SECTIONS 

In order to obtain the average number of collisions and thus the total cross-sections, the 
overall normalization, given by asoft has to be determined. 

As mentioned in the introduction, we are assuming that the entire rise is due to ajet- For 
the proton-proton cross-section, one needs only one further parameter for the non perturba- 
tive region, namely a constant ctq which gives the normalization of the cross-section. As far 
as the proton-antiproton cross-section is concerned, the rapid decrease after the resonances 
is interpreted as dual to the Regge trajectory exchange and it should be described by a 
power • Neglecting the real part of the eikonal, our model is now complete 

and reads as follows : 




(54) 
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FIG. 3: The impact parameter distribution function calculated for the soft gluon summation model, 
using the Qmax values described in the text, for various energy values. Full lines are for the soft 
term. 

with 

x{b, s) = AsNib, + g^l)a,e* (55) 

We also have 



(^Zft = ^0, oZft = ^o(l + — ) (56) 



2 

We find that, in order to properly reproduce the normalization of the cross-section, we need 
a value ctq = 48 m6, in good agreement with the considerations of Sect. II. We now show in 
FIG. the average number of collisions as a function of b, distinguishing between hard 
and soft contributions and using the values of qmax shown in FIG. (0). We only show the 
low energy region, ^/i = 10 100 GeV where the transition between soft and hard processes 
occurs. 

Finally, in FIG. (0) we show the results of our model, putting all the pieces together, for 
the total cross-section for proton-proton and proton-antiproton collisions. We see that the 
model gives an overall satisfactory description of the energy behaviour of available data ^ . 
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FIG. 4: The average number of collisions at -^s = 10, lOOGeV is plotted for the soft gluon 
summation model, using the Qmax values described in the text. The dot-dashed line corresponds 
to the jet contribution at ^/s = 10, 100 GeV. 

IX. ENERGY DEPENDENT < 6^ > in THE BLOCH-NORDSIECK MODEL 

The energy dependent transverse overlap function, Abn, discussed in the previous sec- 
tions, can be used to estimate the energy dependence of the average distance among partons 
in the transverse space during a scattering process, namely 



The energy dependence of the average rms distance between partons so defined, is shown in 



One can see from this figure that for the hard part of the eikonal in the Bloch-Nordsieck 
model (dotted curve), the mean distance between the scattering partons does decrease as 
the energy increases, thereby increasing the shadowing and taming the rise, as opposed to 
the Form Factor model where < 6^ > is a constant. This is then further reflected in a more 
modest high energy rise for the BN model as seen in FIG. (0). It is also pleasing to note 




(57) 
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FIG. 5: Comparison of pp and pp total cross-section data 
from the Bloch-Nordsieck model described in the text. 



Q, Q, Q, Q, Q, B 



with results 



the following self consistency. That is, at low values of ^/s, values of < 6^ > are the same 
for both the soft and hard part of the eikonal in both models, as they must since at low 
energy the transverse overlap function from the BN model is very similar to that from the 
FF model. 

Observations about the need of a shrinkage in the radius of the proton, have been made in 
Ref . 5^, |53| , where multiparticle production in hadron hadron interactions has been studied 
in detail in an eikonal Monte Carlo model. They find that in the hard multi parton model, a 
good fit to the CDF data is obtained if the proton radius is decreased js^ by about a factor 
of 1.7, as compared to the form factor model. Now the observation has been extended to 
photoproduction as well js^. 
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Average rms distance between partons 
in proton-proton scattering 

GRV p,^i„=1.15 GeV, soft from pp, BN 
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FIG. 6: The average distance between partons in the Bloch Nordsieck model (full line) and the Form 
Factor model (dot-dashes) is shown as a function of energy. For the former case, the contribution 
to this distance from the purely soft processes (dashes) with pt < ptmin, and purely hard ones 
(dots) with Pt > Ptmin, is also shown. 

X. CONCLUSION 

In conclusion, we have shown that standard QCD processes such as hard parton parton 
scattering and soft gluon emission from valence quarks can account for two salient features 
of the total proton-proton cross-section, the rise at high energy and the very gentle decrease 
at low energy. An important characteristic of this treatment is that, as the minijet cross- 
section rises with energy, soft gluon emission produces an acoUinearity of the partons and 
reduces the probability of collisions. This affects the cross-sections in two ways : at low 
energy it produces a very soft decrease in afot contributes to the faster decrease in atot, 
at high energy it tames the rise due to a jet- It is then possible to have a very small Ptmin to 
see the onset of the rise around 10 ^ 20 GeV, without encountering too large a cross-section 
when the energy climbs into the TeV range and beyond. We stress that the above behaviour 
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is obtained from leading soft gluon emission from the valence quarks. Subleading emission 
from internal gluon legs is not considered here. It is to be emphasized that singular as 
appears necessary for this purpose thereby implying that confinement plays a crucial role in 
the energy dependence of the total cross section. 

Further input is needed to understand the scale or the normalization, which plays a 
dominant role in the early decrease of the proton-antiproton cross-section. 
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